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Abstract 



The Quark Meson Coupling Model which describes nuclear matter as a collec- 
^ \ tion of non-overlapping MIT bags interacting by the self-consistent exchange 

of scalar and vector mesons is used to study nuclear matter at finite temper- 
Q ■ ature. In its modified version, the density dependence of the bag constant 

^ ■ is introduced by a direct coupling between the bag constant and the scalar 

mean field. In the present work, the coupling of the scalar mean field with 
\ the constituent quarks is considered exactly through the solution of the Dirac 

I equation. Our results show that a phase transition takes place at a critical 

■ temperature around 200 MeV in which the scalar mean field takes a nonzero 

■ value at zero baryon density. Furthermore it is found that the bag constant 
00 ■ decreases significantly when the temperature increases above this critical tem- 

' perature indicating the onset of quark deconfinement. 
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I. INTRODUCTION 



The Quark Meson Coupling Model (QMC), initially proposed by Guichon describes 
nuclear matter as a collection of non-overlapping MIT bags interacting by the self-consistent 
exchange of scalar a and vector uj mesons in the mean field approximation . It is a simple 
extension of the Walecka model except that the meson fields are coupled directly to 
the constituent quarks themselves rather than to the nucleons as in the Walecka model. The 
QMC model thus incorporates explicitly the quark degrees of freedom. This simple model 
was later refined by including the nucleonic fermi motion and the center-of-mass correction 
1^ to the bag energy and applied to a variety of problems . 

In a modification of the original QMC model it has been suggested by Jin and Jennings 
that a reduction of the bag constant in nuclear matter relative to its free-space value 
may be essential for the success of relativistic nuclear phenomenology and that it may play an 
important role in low and medium nuclear physics such as understanding the EMC effect [ITT 



The density-dependence of the bag constant is introduced in the present work by coupling 
it to the scalar meson field as suggested in Ref. [||. It was found by Jin and Jennings that 
when the bag constant is significantly reduced in nuclear medium with respect to its free- 
space value, large cancelling isoscalar Lorentz scalar and vector potentials for the nucleon 
in nuclear matter emerge naturally. Such potentials are comparable to those suggested by 
relativistic nuclear phenomenology and finite density QCD sum rules [|1^ . 



Recently, Panda et al. |TB| have studied nuclear matter at finite temperature using this 
modified version of the QMC model. They determined the scalar mean field by minimiz- 
ing the grand thermodynamical potential with respect to the a field and using the self- 
consistency condition which relates the vector mean field to the baryon density. They found 
that the nucleon properties at finite temperature and/or nonzero baryon density are appre- 
ciably different from their zero temperature vacuum values. 

Our present work is essentially an extension of the work of Panda et al. [|l^ . We intend 
to use the QMC model at finite temperature and to take the medium dependence of the bag 
parameters into account. We however will attempt to solve the self-consistency condition 
for the a field exactly by taking into consideration the full coupling of the scalar mean field 
to the internal quark structure by means of the solution of the point-like Dirac equation 
with the required boundary condition of confinement at the surface of the bag as suggested 
by Refs. . This was not done exactly in the solution of the self-consistency condition for 
the a field by Panda et al. 

The outline of the paper is as follows. In section II, we present the QMC model for nuclear 
matter at finite temperature, together with the details of the self-consistency condition for 
the scalar mean field. In section III, we discuss our results and present our conclusions. 



II. QMC MODEL FOR NUCLEAR MATTER AT FINITE TEMPERATURE 



The QMC model at finite temperature is described in detail in Ref. |jT^. We give here 
the essential equations necessary for the present calculations. The quark field ^^(r, t) inside 
the bag satisfies the Dirac equation 



gluoP ^,(f,t) = 0. 



(2.1) 
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Here the single-particle quark and antiquark energies in units of i? ^ are given as 



where 



(2.2) 



(2.3) 



and m* = — g'^a is the effective quark mass. The boundary condition at the bag surface 
is given by 



9 ' 



(2.4) 



which determines the quark momentum Xnn in the state characterized by specific values of n 
and K. For given values of the bag radius R and the scalar field cr, the quark momentum Xnn 
is obtained from Eq.( p.4|) . The quark chemical potential /ig assuming that there are three 
quarks in the nucleon bag is determined through 



rzg = 3 



3E 
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The total energy from the quark and antiquark is 



(2.5) 



tot 
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R 
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(2.6) 



The bag energy is given by 



bag 



Etot-^ + fR'Bia). 



(2.7) 



where B{a) is the bag parameter. In the simple QMC model, the bag parameter B is taken 
as Bq corresponding to its value for a free nucleon. The medium effects are taken into 
account in the modified QMC model by the following ansatz for the bag parameter |[T3|j3[] 



B = Bq exp 



M, 



(2.8) 



N 



with as an additional parameter. The spurious center-of-mass momentum in the bag is 
subtracted to obtain the effective nucleon mass [HI 



N 



^bag ^ Pcm ^ 



where 



< Pcm > = 



< x^ > 

i?2 



(2.9) 



(2.10) 



and 



3 



UK 
UK 



1 1 



(2.11) 



The bag radius i? is obtained through the minimization of the nucleon mass with respect to 
the bag radius 



0. (2.12) 



dR 

The total energy density at finite temperature T and at finite baryon density p_B is 



(2vr) 



I d'k^/k^Tm/ifB + fs) + + 2^'^'' (2-13) 



where 7 = 4 is the spin-isospin degeneracy factor and fs and f ^ are the Fermi-Dirac 
distribution functions for the baryons and antibaryons 



and 



with e* = \lk? + M^"^ the effective nucleon energy and jj,*^ = jj, — g^uo the effective baryon 
chemical potential. The chemical potential for a given density pb is determined by 

PB = j^Jd'Hh-jB) (2.16) 

where 

00 = ^PB. (2.17) 
The pressure is the negative of the grand thermodynamic potential and is given by 

P = It^, I d'AfB+7B) + Iralco' - l^W- (2.18) 



3(27r)3y e*''" ' ^ ' ^ 2 

The scalar mean field a is determined through the minimization of the thermodynamic 
potential or the maximizing of the pressure §^ = |jl3| . The pressure depends explicitly on 
the scalar mean field a through the last term in Eg. ( ^.1^1 ). It also depends on the nucleon 
effective mass M'^ which in turn also depends on a. If we write the pressure as a function 
of and a [^,0, the extremization of P{M^,a) with respect to the scalar mean field a 
can be written as 

dP _( dP\ (dP\ _ 

da-[dM*J^^, da ^U^Mr ^ ^ ^ 

where 
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(2.20) 



and 



dP 



7 1 



dM: 



3 (27r)3 
7 1 



fB + fl 



3 (27r)3T 
7 1 1 
3 (27r)3T 

duj 



A;2 MJ, 



duj \ 



fB{l-fB) + fBa-fB] 



9u 



dM: 



d'k- 



fB{i-fB)-fBa-f 



Bj 



dM 



(2.21) 



Since the baryon chemical potential and temperature are treated as input parameters, 
the variation of the vector mean field uo with respect to the effective nucleon mass M^ at a 
given value of the baryon density pb reads 



' du \ 





'JbH 


-Jb) 


-fBi^-fB) 


^ + m2 (27r)3 / ^ 


/b(1 


-Jb) 


+ 1b{^-1b) 



(2.22) 



The coupling of the scalar mean field a with the constituent quark in the non-overlapping 
MIT bag through the solution of the point like Dirac equation should be taken into account 
to satisfy the self-consistency condition. This constraint is essential to obtain the correct 
solution of the scalar mean field a. The differentiation of the effective nucleon mass 
with respect to a gives 



dM\ 



N 



da 



TP dEfjgg 



1 1 d<x'^> 

2 B? da 



where 



da 



E 



dE, 



bag dQ'^'^ 



da 



+ 



'dE, 



bag 



da J 



(2.23) 



(2.24) 



and 



d < x"^ > 
d^ 



dni;^ 



' d < 
d^ 



(2.25) 



The depends on x 



Its evaluation can be obtained from the solutions of the point like 
Dirac equation for the constituent quarks which satisfy the required boundary condition on 
the surface of the bag 0,0]. In our case it reads p|,p| 



~d^ 



-glR < ir'\r^ > 



(2.26) 



5 



III. RESULTS AND DISCUSSIONS 



We have studied nuclear matter at finite temperature using the modified quark meson 
couphng model which takes the medium-dependence of the bag into account. We choose a 
direct coupling of the bag constant to the scalar mean field a in the form given in Eq. (|2.8| ). 
The bag parameters are taken as those adopted by Jin and Jennings where B^^ = 188.1 
MeV and Zq = 2.03 are chosen to reproduce the free nucleon mass Mn at its experimental 
value 939 MeV and bag radius Rq = 0.60 fm. The current quark mass is taken equal to 
zero. For = 1, the values of the vector meson coupling and the parameter as fitted from 
the saturation properties of nuclear matter, are given as g'^/2n = 5.24 and g^^ /Atc=3.Q9. 

We first solve Eqs.( |2.16D and (|2.17| ) for given values of temperature and density ps to 
determine the baryon chemical potential fiB- This constraint is given in terms of the effective 
nucleon mass which depends on the bag radius R, the quark chemical potential fig and 
the mean field a. For given values of a and u the bag radius and the quark chemical potential 
fig are obtained using the self-consistency conditions Eq.( |2.12| ) and Eq. , respectively. 
The pressure is evaluated for specific values of temperature and fiB which now becomes an 
input parameter. We then determine the value of a by using the maximization condition 
given in Eq. (|2.19|) . This method is analogous to the method used by Saito and Jin for 
the zero temperature case. It takes into account the coupling of the constituent quarks with 
the scalar mean field in the frame of the solution of the point like Dirac equation exactly. It 
differs from the minimization method used in Ref. |T^. 

Fig.l, displays various pressure isotherms vs. the baryon density. The pressure has the 
usual trend of increasing with temperature and density. It is important to note that the 
pressure attains a nonzero value at zero baryon density above a critical temperature ~ 200 
MeV. This occurs because the scalar mean field a attains a nonzero value at zero baryon 
density Pb = [see discussion concerning Fig. 2 below] as was also observed in the Walecka 
model where it leads to a sharp fall in the effective nucleon mass at Tc. This rapid fall 
of with increasing temperature resembles a phase transition when the system becomes 
a dilute gas of baryons in a sea of mesons and baryon-antibaryon pairs. 

In Fig. 2, we display the scalar mean field cr as a function of the baryon density for various 
temperatures. Fig. 2 (a) indicates that the value of a initially decreases with increasing 
temperature for temperatures less than 200 MeV. However, by the time the temperature 
reaches 150 MeV there are indications of an increase in a at very low baryon densities with 
a nonzero value at pb = 0. For still higher temperatures, as can be seen in Fig.2 (b), the 
situation is more dramatic with the value of a increasing with temperature for all values of 
Pb- This is an indication of a phase transition to a system of baryon-antibaryon pairs at 
very low densities as mentioned above. 

The density and temperature dependence of the baryon effective mass is shown 
in Fig 3. For low baryon density pb, as the temperature is increased, first increases 
slightly and then decreases rapidly for T = 200 MeV for densities less than about 0.2 fm~^. 
This rapid decrease of with increasing temperature resembles a phase transition at a 
high temperature and low density, when the system becomes a dilute gas of baryons in a 
sea of baryon-antibaryon pairs |0].This behavior is consistent with the Walecka model 
and resolves the contradiction that appeared in the earlier calculations. Below the critical 
temperature the effective mass grows with temperature. Above the critical temperature, a 
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increases with temperature thus reducing the nucleon effective mass M^. 

Finally, we display the density dependence of the bag constant for different values of the 
temperature in Fig. 4. The bag constant as shown in Fig. 4 (a), grows with temperature for 
temperatures less than 200 MeV (except at densities smaller than 0.1 fm~^ where B starts 
to decrease for temperatures greater than 150 MeV). However, the situation is completely 
reversed after the phase transition takes place. This is displayed in Fig.4 (b), where the bag 
constant displays a dramatic decrease with temperature for all densities at temperatures 
greater than 200 MeV. This indicates the onset of quark deconfinement above the critical 
temperature: at very high temperature and/or density the hadrons will dissolve into a 
quark-gluon plasma through what is believed to be a first order phase transition. Because 
the QMC model, despite its limitations and shortcomings, uses the quark degrees of freedom 
explicitly, it has made it possible to observe the quark deconfiment phase transition which 
would not have been possible if only the nucleonic degrees of freedom are used. It remains 
to be seen if these calculations can also be carried out with more sophisticated nucleonic 
models and ultimately corroborated by QCD calculations. 
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FIGURES 



FIG. 1. The pressure for nuclear matter as a function of the baryon density ps for various 
values of temperature. 

FIG. 2. The mean scalar field a as a function of the baryon density pB, (a) for different 
temperatures up to T=200 MeV, (b) for different temperatures above T=200 MeV. 

FIG. 3. The effective micleon mass for nuclear matter as a function of the baryon density 
Pb for different values of temperature. 

FIG. 4. The bag constant versus the baryon density pB, (a) for different temperatures up to 
T=200 MeV, (b) for different temperatures above T=200 MeV. 
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Figure 1 




10 



Figure 2a 
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Figure 3 
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Figure 4a 
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